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Abstract :  We  analyze  the  robustness  of  various  standard  finite  element 
schemes  for  the  Reissner-Mindlin  plate  and  obtain  asymptotic  convergence 
estimates  that  are  uniform  in  terms  of  the  thickness  d.  We  identify  h 
version  schemes  that  show  locking,  l.e.  for  which  the  asymptotic  convergence 
rate  deteriorates  as  d  ->  0  and  also  show  that  the  p  version  is  free  of 
locking.  In  order  to  isolate  locking  effects  from  boundary  layer  effects 
(which  also  arise  as  d  -»  0),  our  analysis  is  carried  out  for  the  periodic 
case,  which  is  free  of  boundary  layers. 
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1. 


INTRODUCTION 


It  is  well-known  that  the  numerical  approximation  of  the  Reissner-Mindlln 
(RM)  plate  model  by  certain  finite  element  schemes  deteriorates  when  the 
thickness  d  of  the  plate  is  close  to  zero.  This  occurs  because  of  two 
phenomena,  the  existence  of  boundary  layers  and  the  presence  of  locking. 

Boundary  layers  arise  as  components  of  the  exact  solution  of  the  RM  plate 
model  and  (to  some  extent)  model  the  boundary  layers  present  in  the  exact 
solution  of  the  corresponding  3-d  problem.  In  a  series  of  papers  (see  [1]  for 
example),  Arnold  and  Falk  have  characterized  in  detail  various  aspects  of  the 
boundary  layer  for  the  exact  RM  solution,  for  different  boundary  conditions. 

In  particular,  in  [1]  it  is  shown  that  the  layers  are  strongest  for  soft 
simply  supported  and  free  plates  and  weakest  for  soft  clamped  conditions. 

One  effect  of  the  presence  of  these  boundary  layers  is  to  weaken  the  a  priori 
regularity  of  the  solution.  Since  the  singular  behavior  occurs  only  near  the 
boundary,  an  effective  strategy  to  overcome  any  consequent  deterioration  of 
numerical  schemes  is  to  locally  refine  the  mesh  (usually  with  a  low  order 
scheme).  See,  for  Instance  [2],  where  this  has  been  done  in  the  context  of 
the  h-p  version. 

The  second  phenomenon  mentioned  above,  numerical  locking,  also  occurs 
when  the  thickness  is  close  to  zero,  but  for  a  different  reason,  it  is 
well-known  that  as  the  thickness  d  -»  0,  the  RM  solution  tends  to  the 
Klrchoff  plate  solution.  In  the  limiting  case,  the  exact  solution  will 
therefore  satisfy  Kirchoff’ a  hypothesis.  The  problem  of  locking  occurs 
because  as  d  0,  the  finite  element  solution  is  also  forced  to  satisfy  this 
hypothesis.  Consequently,  the  number  of  conforming  trial  functions  (which 
satisfy  Klrchoff’ s  hypothesis)  can  get  severely  restricted,  resulting  in  a 
degradation  of  the  approximation  properties  of  the  trial  space. 

One  method  to  avoid  locking  is  to  construct  a  finite  element  space  which 
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possesses  optimal  approximation  properties  even  when  restricted  by  Klrchoff*s 
constraint.  This  leads  to  a  robust  standard  finite  element  method,  i.e.  one 
whose  performance  is  not  sensitive  to  the  thickness.  An  alternate  strategy  is 
to  use  a  mixed  method,  which  has  the  effect  of  enforcing  the  constraint  in  a 
weaker  sense,  thereby  restricting  the  finite  element  space  less.  Such  methods 
have  been  developed  for  various  contexts  where  locking  occurs  —  see,  for 
example,  [9]  (particularly  Chapter  7  for  the  R-M  plate)  and  the  references 
therein.  In  this  paper,  we  will  concentrate  exclusively  on  the  first  strategy, 
which  can  be  used  without  reformulating  the  usual  variational  form.  As  a 
consequence,  our  results  will  be  Immediately  applicable  to  several  finite 
elements  used  in  the  context  of  various  commercial  codes. 

The  problem  of  locking  is  quite  different  in  terms  of  origin  and 
numerical  treatment  from  that  of  boundary  layers.  Consequently,  it  is  more 
instructive  to  analyze  these  problems  separately.  In  this  paper,  we  will  be 
Interested  only  in  the  problem  of  locking.  To  isolate  this  phenomenon  and 
divorce  it  from  the  effect  of  boundary  layers,  we  will  be  considering  the  case 
of  periodic  boundary  conditions,  which  we  will  choose  so  that  the  solution  is 
smooth.  (This  models  the  situation  in  the  interior,  in  the  case  boundary 
layers  are  present. )  For  this  RM  model  problem  we  will  characterize  the 
locking  and  robustness  properties  of  various  finite  element  schemes,  using  the 
general  theory  of  locking  developed  by  us  in  [4].  A  key  condition  from  that 
work,  the  so-called  "condition  (a)",  will  be  shown  to  be  satisfied  in  this 
case,  thereby  reducing  the  question  of  locking  to  one  of  approximability 
alon«*  A  similar  technique  was  used  by  us  in  [5]  to  analyze  Poisson  ratio 
locking  which  occurs  in  elasticity  problems  when  the  Poisson  ratio  is  close  to 
1/2. 
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2.  THE  MODEL  PROBLEM  AND  ITS  REGULARITY 

We  consider  as  our  domain  the  2-d  midsection  Q  3  (-n.it)  of  a  square 
isotropic  plate  with  the  plate  occupying  the  region  Q  x  [-d,d].  On  Q  we 
consider  the  RM  plate  model  for  u^  = 

(2.1)  Cdud  »  -d2  5  {(l-v)  +  (1+w)  77*^d>  ‘  A(Vud-^d)  *  0 

<2. 2)  -Ad‘27*(7«d-^d)  *  g 

which  gives  the  bending  of  the  plate  in  equilibrium.  We  assume  the  periodic 
boundary  conditions, 

(2.3a)  <*>d ( x , n )  *  w^x.-ir),  ud(rt,y)  =  ud(-n,y),  |xj,|y|  s  n 

(2.3b)  ^d(x,H)  =  ^d(x,-it),  ^d(H,y)  =  ?d(-H,y),  |x|,  I  y  I  Sir. 

Here,  $d  gives  the  rotation  of  fibers  normal  to  the  midplane,  u>d  measures 

E 

the  transverse  displacement  of  the  midplane,  and  D  *  - j-  is  the 

12(l-v  ) 

3  kE 

flexural  rigidity  scaled  by  d  .  Also,  A  »  2(l*v)  wlth  K  >  °>  E  and  0  s 
v  <  0.5  being  the  shear-correction  factor.  Young’s  modulus  and  Poisson  ratio, 
respectively.  Moreover,  g  is  the  scaled  loading  function,  l.e.  the 

3 

transverse  load  density  per  unit  area  divided  by  d  .  (We  assume, 

3 

essentially,  that  the  loading  function  for  a  plate  of  thickness  2d  is  gd  , 
where  g  is  independent  of  d. ) 

If  we  allow  d  to  tend  to  zero  in  (2.1),  we  formally  obtain  Kirchoff’s 
hypothesis, 

(2.4)  CQu0  :=  $Q  ~  grad  »Q  =  0 

The  limiting  solution  uQ  =  (^0,uQ)  satisfies  (2.4)  together  with 

(2.5)  DA2o>0  =  g 
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and  the  periodic  boundary  conditions  (2.3)  (see  e.g.  [31).  This,  of  course, 
is  the  classical  blharmonic  equation  of  plate  bending,  which  (unlike  the  RM 
model)  is  independent  of  d  (and  is  sometimes  used  to  model  the  actual  3-d 
plate). 

n  s 

For  any  domain  R  c  R  .  we  will  denote  by  H  (R)  the  usual  Sobolev 

space  with  s  (=  0,1,2,...)  square  lntegrable  derivatives.  For  our  domain 

Q,  Hs  (£3)  will  denote  the  corresponding  space  of  functions  with  s 
per 

periodic  derivatives,  periodic  in  both  x  and  y.  Using  the  method  of  real 

interpolation  [71,  the  above  spaces  may  be  defined  for  all  real  s.  We  will 

use  j»|  I*IS  to  denote  the  seminorm  and  norm  respectively,  in  both  the 

(s) 

periodic  and  non-periodic  case.  We  will  also  use  (^^(£3),  which  will 

denote  the  space  of  functions  with  s  periodic  continuous  derivatives. 

Any  u(x)  »  u(x, ,x„)  in  HS  (£3)  can  be  expanded  as  a  Fourier  series, 
1  JL  per 


(2.6) 

where 


u(x) 


£  ukeik*x,  k  =  (krk2).  x  =  (xrx2) 
keZ2 


k 


u 


with  <•,*>  denoting  the  usual  L2(£3)  inner  product.  Then  we  have  the 
following  equivalences: 


(2.7) 


*  l  |k|2Vl2.  M2 

keZ2 


l  (l+|k|2)s|uk!2  , 
keZ2 


which  hold  for  all  real  s  a  o.  We  may  also  use  (2.7)  to  define  the  norms 

g 

|ujs,  s  <  0  and  the  corresponding  spaces  ((^^(£3)  which  are  the  completions 

under  these  norms  of  (£3). 

per 

We  now  cast  (2.1)-(2.3)  into  the  following  weak  form.  Find,  for 
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d  e  (0,1],  u .  *  (^.,w.)  e  V  =  [H1  (fl)]3  such  that  for  all  v  *  (t, €)  e  V, 

a  aa  per 


(2.8) 


ad(Vv)  =  b(ud*v)  +  ^2  <coud-cov>  "  F(V) 

a 


Here,  CQ  is  as  defined  in  (2.4), 

b(ud,v)  =  \  JJ  <<l-v)V?d*vt  +  (l-*-y)(v4d)(v4)>dx1dx2 

and 

F(v)  =  <g,£>  =  JJ  dx1dx2 

with  the  last  integral  being  understood  as  the  pairing  of  H1  r(£l)  and 

H  *  (£2)  if  g  e  H  *  (Q).  (For  existence  and  uniqueness  for  this  problem,  see 
per  per 

Theorem  2.1  ahead.)  We  will  denote  the  problem  (2.8)  by  Prf  and  assume, 
without  loss  of  generality,  that  it  is  equivalent  to  (2.1)-(2.3). 

Let  us  define  the  energy  norm  corresponding  to  P^  by 

M^d  =  ad{u’u)’  ueV‘ 

Also,  for  u  e  H*  (0)  x  H*  (£1)  x  H*  (fi)  =  H^fQ),  let 

per  per  per  per 

Mk.t  -  -  \*i\l  *  1*2*1  *  l“lf  • 

where  the  norm  in  u  is  understood  to  be  modulo  constants.  Then  it  may  be 

shown  that 


(2.9) 


AiNi,i  s  ME,d  s  A2d’  Hi,i 


where  and  are  constants  Independent  of  d.  It  is  readily  observed 

that  for  d  bounded  away  from  0,  the  two  norms  are  equivalent. 

Jr+1  Jr 

Let  us  now  define,  for  1c  a  1,  0  s  d  <  1,  the  spaces  H.  .  c  m  ’  (Q) 

it,  d  per 


given  by 


(2.10) 


Vd  ’  cd“  =  °>- 
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For  any  normed  linear  space  H,  we  will  denote  the  ball  of  radius  B  >  0  by 

H®  =  (u  €  H,  |u|H  s  Bh 

We  will  use  the  notation 

4  3  -  (u  *  H^*’k(a),  Mk,lk  »B) 

and 

<<«  -  <  "  "k.r 

We  note  that  (2.10)  may  equivalently  be  characterized  as 


0} 


since  using  (2.1),  it  is  easily  seen  that  %  e  Hlt+1(£2)  =*  u  e  H*  (£1). 

per  per 

Similarly,  in  the  definition  of  we  may  replace  Suijc+i  ^  by 

We  note  also  that  for  u  €  d>  we  have,  using  (2.1), 


k+1 


(2.11)  lC0Ulk-l  3 

for  K  a  constant  Independent  of  u,  d. 

Let  us  look  more  closely  at  the  limiting  sets  H^  Q  and  Q.  Since 

for  these  CLu  =  0,  we  see  that  |  €  H*+*(Q)  implies  that  Vu  =  ^  e  Hlc+1(n), 
u  per  per 

so  that  w  e  H^+^(Q) .  Hence,  we  see  that  in  this  limiting  case,  the 
per 

regularity  of  w  is  increased  by  two  derivatives  and  we  have  the  equivalent 
character i zat ion , 


Vo  3  <i«*&J’k+2(n).  c0u  =  0) 

’’k.O  (u  *  Vo’  luL+l,k+2  3  B> 

The  choice  of  the  above  is  motivated  by  the  following  theorem,  which 
gives  an  a  priori  estimate  in  these  weighted  norms. 
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s-2 

Theorem  2. 1 .  Let  g  €  H  (Q),  s  £  1,  satisfy  the  compatibility 

pe  r 

condition 

(2.12)  <g,l>  *  0. 

Then  there  exists  a  unique  sequence  of  solutions  (u^)  *  ^d,Wd^  ^wd  unique 
up  to  a  constant)  to  (2- 1 )— (2. 3)  for  d  e  (0,1)  and  (2.3)~(2.5)  for  d  *  0, 
such  that  ud  €  d»  where  B  =  C|8|s  with  C  a  constant  independent  of 
g  and  d. 

Proof.  Suppose  g,  0^,  ^2  an<^  wd  are  rePresen^®^  terras  of  their 
respective  Fourier  series,  as  in  (2.6).  Then  (2.1)-(2.2)  may  be  written  as 
(d  >  0) 


^  (  2k2 +  ( 1 -v ) k2 } +Xd_2  gd+tOk^ 

£  eik’X  5  (1+V)k1k2  5{2k2+(l-w)k2>+Xd2 

keZ2 

iAd”2k1  iAd2k2 


where  Zq  *  ZMO}.  Here,  we  have  used  the  fact  that  g°  *  0,  due  to  (2.12). 
Similarly,  the  solution  of  (2.3)-(2.5)  is  given  by 
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(2.14) 


k 


ik*z  g 


Dlkl 


with  =  7uq.  Note  that  In  (2. 13 ) — ( 2 . 14) ,  Is  unique  up  to  a  constant. 


From  (2.13),  It  is  easily  seen  that  for  d  >  0, 


l"dl,M.s  *  Cl*IS-2 

proving  that  a  ^  .  The  case  d  =  0  follows  similarly  from  (2.14).  □ 

We  now  prove  the  following  theorem,  which  gives  us  the  co-called 

“condition  (a)"  which  is  central  to  the  locking  theory  developed  in  the  next 

section.  Essentially,  this  condition  establishes  the  rate  at  which  the 

solution  of  (2. 1 ) — ( 2 . 3 )  tends  to  the  solution  of  (2.3)-(2.5)  in  the  HS+1‘S((1) 

per 

norm  when  d  is  small. 


Theorem  2.2.  For  any  u,  =  (^..w.)  e  H5  .,  there  exists  u„  =  ($n,un)  € 

d  d  d  s.d  u  uu 

HcB_  such  that 
s,0 


l“d-“oU  3 

where  c,K  are  constants  independent  of  u^,  d  and  B. 


Proof.  Given  u.  =  (^.,w.)  e  .,  let  us  define  g  (possibly 
d  d  d  s,d 

depending  on  d)  by  (2.2).  Then  by  (2.2), 

Ws-2  S  Kd'2|C0“dls-l 

s  KB 

where  we  have  used  (2.11). 

Now  may  be  expressed  in  terms  of  g  by  (2.13).  Also,  let 

be  given  by  (2.14)  and  define  =  VwQ.  Then  it  is  easily  seen  that 
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h  *  t 


0* 


■A  £g 

and  that  (?q,«q)  €  Hg  Q,  where  c  is  independent  of  ud,  d  and  B. 
tuting  *  ^w0  in  ^2-^*  we  have 

VC«d-«0)  =  -A_1d2  5  {(i-v)A?d  +  (l+u)  7V*^d> 


Substl- 


froa  which  it  follows  that 


l“d '  "oU  3  Kd2l?dU*i  3  m2b- 


The  above  theorem  shows  that  for  d  small,  functions  in  H®  ,  are  close 

s,d 

to  functions  in  H®  _.  In  this  connection,  we  will  also  need  the  following 

S  |  U 

result. 


Theorem  2.3.  Given  uQ  €  q,  s  a  1,  there  exists  a 
constant  C  independent  of  uQ,d  such  that  for  any  d  e  (0,11,  there 
is  a  ud  €  Hg  d  satisfying 


(2.15) 
and 

(2.16) 


Kkd  *  l“dU.S  3  Cluols*l.s*2 


•ud  '  uoll,l  "*  0  33  d  *  °- 


Proof.  First,  using  (2.11),  it  is  easily  seen  that  for  s  a  1, 


(2.17) 


ludlE,d  *  CIud^s+l,! 


Next,  let  u_  *  (grad  w_,«_)  e  H®  ...  Then,  since  u>n  «  HS+2(Q),  we  obtain  the 
u  oo  s,o  o  per 

decomposition  (wQ  is  defined  modulo  constants) 


“0  ■  I 

k*Z^ 


e 


uk.  j;  (|k|2)s*2|0ki2  <  .. 

k3*0 
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k  k  4 

Hence,  defining  g  =  Dw  |k|  ,  we  see  that 

(2.18)  8«H^(0),  <g,l>  .  0,  |g|s.2  ,[Uois.2  *  C|u0[st]  ^  . 

Using  Theorem  2.1  together  with  C2. 17)— <2. 18)  allows  us  to  construct  ud 
satifying  (2.15).  Also,  the  argument  of  Theorem  2.2  shows  that  (2.16)  holds. □ 


3.  LOCKING  AND  ROBUSTNESS 

Suppose  now  that  we  are  interested  in  approximating  (2.8).  We  assume 

that  we  are  given  a  sequence  <V^>  of  finite-dimensional  subspaces  of  V  = 
1  3 

[H  (fl) ]  (N  denoting  the  dimension,  N  €  AO.  Then  we  can  define  the 
per 

sequence  of  finite  element  solutions  u^  €  v"  by 

(3.1)  ad(ud,v)  =  ad(ud,v)  V  v  e  V1**. 


The  sequence  <V*>  thus  defines  an  extension  procedure  ?,  i.e.  a  rule  by 
which  we  can  increase  the  dimension  N  with  the  idea  of  increasing  accuracy. 
(3.1)  immediately  gives 


(3.2) 


iud  -  udlE,d  a  ln£  lud  -  vkd 

veV 


As  shown  in  Theorem  2.1,  depending  upon  the  regularity  of  the  data  g, 

Jc+1  k 

the  exact  solutions  of  our  problem  will  belong  to  the  sets  d  c  Hper'  fQ), 
k  a:  1,  introduced  in  Section  2.  We  will  assume  that  the  sequence  (V^>  is 


Fg-admissible,  i.e.  it  leads  to  a  certain  fixed  rate  Fq(N)  of  convergence 
when  functions  in  HkVl.k 

(Q)  are  approximated,  in  the  following  sense: 

per 


(3.3)  A1F0(N)  *  sup  inf  B u_v R i  i  s  A2F0(N)- 

ueH^  ve/ 

Hero,  Fq ( N )  -»  0  as  N-»  «  and  A^,  depend  upon  B  but  are  independent 

of  N.  Moreover,  we  assume  that  there  exists  d^  e  (0,1)  such  that  for 
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dg  ^  d  S  1 , 

(3.4)  A^d^TgCN)  a  sup  inf  lu-v^  ^  s  A2(dQ)F0(N). 

ueH^  d  veV* 


(Note  that  the  lower  bound  in  (3.3)  follows  from  the  one  in  (3.4)  while  the 
upper  bound  In  (3.4)  follows  from  the  one  in  (3.3)). 

Using  (2.9),  (3.2)  and  (3.4),  we  then  obtain  the  following  estimate  which 
holds  uniformly  for  all  d^  s  d  s  1, 


(3.5) 


where 

(3.6) 


N, 


WF0(N)  *  SUP  Ed(ud*V  S  *2(d0)FO(N) 

V<d 


Ed(v)  -  H,_,  or  |v|E  d 


Whether  or  not  Xj,  X2  are  bounded  as  d^  -»  0  will  depend  on  the  extension 
procedure  being  used. 

A  procedure  9  for  which  (3.5)  holds  uniformly  for  ail  0  <  d  s  l  will 
be  called  free  from  locking  for  the  sets  d  with  respect  to  the  Ed 
measure.  A  more  precise  definition,  adapted  from  the  general  treatment  of 
locking  in  14],  is  given  below. 

Let  L(d,N),  the  locking  ratio  corresponding  to  d  e  (0,1],  N  e 
with  respect  to  the  spaces  d  c  H^^,,C(n)  and  error  measures  (Ed>  (as  in 
(3.6))  for  the  problems  (3.1),  be  defined  by 

L(d,N)  =  sup  Ed(ud-ud)(F0(N))'1. 


Then  we  make  the  following  definitions. 


Definition  3. 1 .  The  extension  procedure  f  is  free  from  locking  for  the 

family  of  problems  (3.1),  d  €  (0,1]  with  respect  to  the  solution  sets 

H.  „  c  Hk*1,lt(D)  and  error  measures  E.  if  and  only  if 
K,d  per  4 


12 


■  M  <  to. 


lim  sup  sup  L(d,N) 

l  defo.ll  J 


N-*» 


3F  shows  locking  of  order  f(N)  if  and  only  if 


0  <  11b  sup  jsup  L(d,N)  j  *  K 

N-xa  ^  j 


<  a 


where  f(N)  -*  «  as  N  -»  ®.  It  shows  locking  of  gjt  least  (respectively  at 
most)  order  f(N)  if  K  >  0  (respectively  K  <  •). 

Definition  3.2.  The  extension  procedure  f  is  robust  for  the  family  of 

problems  (3.1),  d  e  (0,1]  with  respect  to  the  solution  sets  H,  .  c  H*+1,lc(Q) 

k,d  per 

and  error  measures  if  and  only  if 


Lim  sup 
N-*»  d 


sup 

V^.d 


0. 


It  is  robust  with  uniform  order  g(N)  if  and  only  if 


sup 

d 


sup 

v£d 


Ed(Vud> 


S  g(N) 


where  g(N)  0  as  N  ->  ». 


Remark.  In  Sections  4,  5  we  will  frequently  use  the  form  g(N)  *  N-r  to 
characterize  the  robustness  order  r.  If  r  »  0,  then  convergence  will  not  be 
guaranteed. 


Definitions  3.1,  3.2  are  related  by  the  following  theorem,  from  [4]. 

Theorem  3.1.  9  is  free  from  locking  if  and  only  if  it  is  robust  with 

uniform  order  Fq(N).  Moreover,  suppose  f(N)  is  such  that 

f(N)FQ(N)  -  g(N)  0  as  N  ->  «. 


13 


Then  ?  shows  locking  of  order  f(N)  if  and  only  if  it  is  robust  with 
maximum  uniform  order  g(N). 


It  is  easily  seen  that  ?  is  non-robust  if  and  only  if  it  shows  locking 
of  order  (FQtN))  1. 

Let  us  briefly  explain  the  above  ideas.  We  are  assuming  that  our  exact 

solution  has  a  certain  regularity  (i.e.  it  is  in  H*+1,lc(Q)).  Our  extension 

per 

procedure  ?  has  associated  with  it  a  rate  of  best  approximation  Fq(N), 
which  gives  the  best  approximation  that  we  could  expect  to  achieve  (Eqn. 

Jr+1  lr 

(3.3)),  for  the  most  unfavorable  exact  solution  in  Ii  ’  (Q).  Using  the 

per 

finite  element  method  gives  this  rate  for  d  e  d^  (Eqn.  (3.4)),  but  does  not 
necessarily  give  this  rate  (uniformly)  as  d  approaches  0.  The  locking 
ratio  compares  the  accuracy  actually  achieved  by  the  finite  element  method 
(for  the  least  favorable  exact  solution  in  ^)  to  the  best  accuracy 
possible  (i.e.  to  FQ(N)).  If  the  achieved  accuracy  is  asymptotically  the  same 
as  the  rate  FQ(N)  in  (3.3),  (3.4)  uniformly  for  all  d  e  (0,1],  then  we  say 
that  the  extension  procedure  i.  is  free  of  locking  for  all  ud  e  d> 

d  e  (0,1].  If  the  achieved  accuracy  is  asymptotically  not  the  same,  then  the 
robustness  g(N)  gives  the  best  rate  of  convergence  that  can  be  achieved 
Independent  of  the  parameter  d.  In  this  case,  f(N)  *  g(N)(FQ(N))  * 
characterizes  the  asymptotic  strength  of  the  locking. 


In  [4],  we  have  formulated  a  useful  condition,  called  condition  (a), 
under  which  the  question  of  locking  reduces  to  one  of  approxlmabllity  alone. 
For  this  condition  to  hold,  we  must  first  be  given  a  sequence  of  solution 
spaces  d  and  a  limit  space  Q  such  that  Theorem  2.3  holds.  Then 

there  exists  a  uQ  e 

(for  some  c  Independent  of  u^,  d,B;  u^  depending  on  ud)  such  that 


condition  (a)  may  be  stated  as:  Given  u.  e 


< 


d’ 


(3.7) 


Jud  *  unl 


0»k+l  ,k 


KdB 


14 


with  K  a  constant  independent  of  B,  d  and  u^.  This  condition  therefore 

characterizes  the  distance  of  solutions  u^  to  functions  u^  satisfying 

(2.4),  as  d  -»  0.  We  have  shown  in  [4]  that  if  the  “remainder’1  u,  -  u-  is 

d  0 

small  in  the  sense  of  (3.7),  then  we  need  only  consider  the  approximation  of 
functions  u^  in  the  limit  space  ^  to  answer  questions  about  locking  and 
robustness. 

Our  choice  of  periodic  boundary  conditions  for  the  plate  problem  is 

motivated  primarily  by  the  fact  that  condition  (a)  is  satisfied.  To  prove 

(3.7)  for  our  problem,  we  simply  choose  u^  (for  given  u^)  as  in  Theorem 

2 

2.2.  Then  we  get  (3.7);  in  fact,  we  get  a  power  of  d  (Instead  of  just  d, 
as  needed).  As  a  result.  Theorem  2.4  from  [4]  will  hold  for  our  problem. 

Tills  theorem  states  that  locking  and  robustness  rates  are  the  same  no  matter 
which  error  measure  in  (3.6)  is  used.  It  is  stated  below. 


Theorem  3.2.  Consider  the  family  of  problems  (3.1),  d  e  (0,1]  with  the 

solution  sets  d  c  Hj^*,lc(fl),  k  s  1.  Then  the  extension  procedure  9  is 

free  from  locking  with  respect  to  the  V  =  H1,*(Q)  norm  if  and  only  if  it  is 

per 

free  with  respect  to  the  energy  norm.  It  shows  locking  of  order  f(N)  in  the 
V  norm  if  and  only  if  it  shows  locking  of  order  f(N)  in  the  energy  norm. 


We  will  now  only  refer  to  the  locking  of  f,  without  specifying  an  error 
measure.  The  following  theorem  reduces  the  question  of  locking  to  one  of 
approximability  alone. 


Theorem  3. 3.  Consider  the  problems  (3.1)  with  solution  sets  d, 

2 

k  fc  1.  Let  V**  *  Y**  x  ^  where  Y**  c  [h1  (fl)l  ,  H1  (£1)  and  define 

L  Per  J  per 

W^cZ^  byWN={weZN,  grad  u  €  Y^>.  Then  the  extension  procedure 
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£  s  {v^>  is  robust  with  uniform  order  max{FQ(N) ,g(N) }  where  g(N)  is 
given  by 

(3.8)  g(N)  *  sup  inf  (<•>  -  z|_  . 

wmH**2'8  zeW* 
per 

Also,  with  Fq(N)  as  in  (3.3),  5  is  free  from  locking  if  and  only  if 

(3.9)  g(N)  s  CFq(N). 

It  shows  locking  of  order  f(N)  if  and  only  if 

(3.10)  C  FQ(N)f (N)  s  g(N)  s  C2FQ(N)f(N). 


Proof.  Let  us  define 


(3.11) 


g(N) 


sup  inf  |u  -  v|  . 

*•> u  n 


cov«o 


Since  condition  (a)  is  satisfied,  by  Theorem  2.2(B)  of  [4],  9  is  robust  with 
uniform  order  max{FQ(N) ,g(N) }.  The  argument  is  as  follows.  For  ud  «  d, 
we  may  find  (by  condition  (a))  a  uQ  €  Q  such  that  (3.7)  holds.  Then  we 
have 

l“d  '  "ali.i  a  l"d  '  “dkd  a  lnf„  l“d  ‘  vlE.d 

vev 


{|u° '  Vl'E’d  *  |lUd*uo) ' 


1  2 
C0V1*° 


.-1 


s  g(N)  ♦  d  inf  |(ud-uQ)  -  v2|1(1 


vX 
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using  (3.3)  and  (3.7)  (the  latter  giving  (u^-Uq)  e  H^0).  This  proves  the 
robustness  order.  Next,  using  the  fact  that  Theorem  2.3  and  condition  (a) 
hold,  we  see  by  Theorem  2.2(B)  of  [4],  that  (3.9)  and  (3.10)  will  hold  with 
g(N)  replaced  by  g(N). 

To  show  g(N)  and  g(N)  are  equivalent,  consider  a  v  =  ($,z)  e  such 
that  CQv  *  0.  Then 

CQv  =  0  *  zeW^,  $  «  grad  z. 

Hence,  for  u  =  (^,«)  e  g,  we  have 


g(N) 


sup 

u<o 


inf 

veV* 


sup 

««l£*?,B(Q) 

per 


inf  { Jgrad  u  -  grad  zl  +  \u  -zL> 

W 


from  which  it  is  easily  seen  that  g(N)  and  g(N)  are  equivalent.  The 
theorem  follows. 

Let  us  now  define  the  subspaces  'Z'  =  /  x  ZN.  We  will  consider  both 
triangular  and  rectangular  meshes.  Our  results  will  be  for  the  uniform 
triangular  and  rectangular  meshes  shown  for  the  case  of  our  square 

domain  Q  in  Figure  3. 1 . 


Fig.  3.1.  Uniform  meshes  M^,  H^. 


Let  S  be  a  triangle  or  parallelogram.  Then  we  define 


*p<s>  - 

polynomials  on 

S 

of  total  decree  a  D 

?>2(S)  * 

p 

polynomials  on 

S 

of  decree  s  d  in  each  variable 

^(S)  = 

p 

^(S)  .  {xPx2. 

Xj> 

the  serendipity  elements 

(*>2(S), 

p 

?>3(S)  are  defined 
P 

only  for  S  parallelogram). 

For  any  mesh  we  now  define  for  1  =  1,2,3 


K'.-i  *  ’,p:-i(Mh)  =  {v  e  l2(o)-  'ij, 6 


where  n"  are  the  elements  of 


For  k  >  0,  we  define 


^•5  =  n  C^(Q) 

p»  *»  p» «  P)  1  per 


n  v  **> 

where  r’.  and  r\  are  defined  for  rectangular  meshes  only.  In  the 
p , K  P,K 

succeeding  sections,  we  will  consider  the  locking  effects  of  the  spaces 

Vi,h  defined  by  taking  Y*1  *  Y1,h  =  and  Z1*  =  Z*’h  *  **•!?.  Note 

p,q  p  p,0  q  q,(J 

that  for  the  space  N  =  0(h  2p2). 

p,  u 


4.  THE  h  VERSION 

In  this  section,  we  consider  the  robustness  of  the  finite  element  spaces 

V*  -  V*'*1  when  p,  q  are  held  fixed  and  h  Is  decreased  to  attain  accuracy, 
p.q 

Let  us  first  estimate  Fq(N)  In  (3.3),  (3.4). 

Lemma  4. 1 .  Consider  an  h  version  seauence  of  spaces  ,v">.  {vj;*} 

2 

defined  by  v1’*1  ■  Y*'*1  x  z*’*1  *  on  an  appropriate  triangular 

p.q  p  q  L  p-°J  q-° 

or  parallelogram  quasiuniform  mesh»  where  p,q  m  1  are  fixed  and  h  varies. 
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Then  for  any  d^  >  0,  (3.4)  is  satisfied  with 

(4.1)  Fq(N)  »  CN-nin(p,q,k-l)/2 

with  C  independent  of  N  but  depending  on  p,q,k.  Moreover,  (3.3)  is  also 
satisfied  with  the  above  Fq. 

Proof.  We  first  note  that  for  l  >  1,  r  a  1,  with  3  *  mln(r,£-l )/2, 

(4.2)  C  »‘?  s  sup  inf  |z  -  v|  s  C j'0  - 

(4.2)  has  been  shown  in  Theorem  4.1  of  [5]  (for  example)  for  the  non-periodic 
case,  from  which  the  periodic  case  is  easily  deduced.  Using  (4.2),  it  is 
easily  seen  that  the  upper  estimate  in  (3.3)  (and  hence  (3.4))  holds  with  F^ 
given  by  (4.1).  To  establish  the  lower  estimate  in  (3.4)  (and  hence  (3.3)), 
we  note  that  given  a  periodic  «d  (or  conversely  a  periodic  ^),  we  can  find 
a  corresponding  periodic  function  %  (or  a  periodic  a>d)  satisfying  (2.1), 

such  that  ud  ■  ^d*ud^  ^or  **d  *  ^d’^d^  lles  ln  ^  d  ^for  aPProPriate 

k).  Hence,  (4.2)  can  be  used  to  establish  the  required  lower  bound.  a 

We  may  now  analyze  the  locking  and  robustness  properties  of  various 

types  of  elements  by  calculating  g(N)  defined  by  (3.8)  and  comparing  it 

with  FQ(N)  given  by  (4.1).  We  first  consider  triangular  elements. 

Suppose  that  V**  *  V*’**  (as  in  Lemma  4.1)  is  defined  on  mesh  Then 

P.q  1 


we  see 

that 

(4.3) 

w" 

s 

q.l 

q  * 

P 

= 

P+1 . 1 

q  a 

P  *  1. 

Let  us 

def lne 

r(q) 

*  max(2,q) 

for  1 

s  q  s  4 
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=  q  +  1  for  q  a  5. 

Then  the  following  lemma  follows  from  Theorem  5.1  and  Lemma  5.2  of  [5], 
where  It  has  been  established  for  the  non-periodic  case. 


Lemma  4.2.  Let  be  the  uniform  triangular  mesh  ot  Figure  3.1. 

for  w  e  H***(Q),  k  a  0.  q  a  1, 
per 


Then 


inf  [<,  -  z|,  *  Ch,ln(k’r(<1,-2)M 


zeP 


l.h 

q,i 


k+2. 


,(«#) , 


Moreover,  there  exists  a  function  Q  e  C  (Q)  satisfying 

per 


inf  |Q  -  z|2  a  Ch 


r(q)-2 


zeP 


l.h 

q.i 


We  can  now  prove  the  following  theorem. 


Theorem  4. 1.  Let  the  extension  procedure  ?  consist  of  the  h  version  on 

2 

the  uniform  triangular  mesh  M^,  with  spaces  <*>  -  {vp:S}  -  {K:S]  x 

Let  the  solution  sets  be  Then  for  the  p  and  q  shown  below,  ?  is 

robust  with  uniform  order  N  r  when  k  a  2r+l  and  shows  locking  of  order  N 
when  k  a  p  +  1. 


Robustness 

Locking  order 

Degree 

Degree 

order  *  0(N  r) 

=  O(N^) 

P 

q 

r 

l 

1 

0 

1/2 

OSB 

q  *  p 

(p-2)/2 

1 

■HI 

q  a  p  +  l 

(p-1 )/2 

1/2 

q  “  p 

(p-D/2 

1/2 

■  1 

q  a  p  +  l 

p/2 

0 

(For  p  >  q,  the  same  results  as  for  the  case  p  *  q  hold. ) 
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Remark.  The  above  theorem  shows  that  with  the  customary  choice  p  =  q, 
convergence  is  not  guaranteed  for  p  =  1  or  2  (i.e.  r,  the  robustness  order 
is  0)  and  is  only  guaranteed  (with  reduced  order)  if  p  £  3.  In  fact, 

locking  cannot  be  avoided  whenever  we  take  p  «  q.  It  is  0(N)  for 

l/?  1  /? 

p  *  2,3,4  and  0(N  )  for  p  £  5  (for  p  =  1,  it  is  technically  0(N  ) 

-1/2 

as  well,  since  the  maximum  possible  rate  in  this  case  is  0(N  )).  To  avoid 

locking,  p  has  to  be  taken  to  be  5  or  larger  with  q  being  chosen  to  be 
p  +  1.  (Note  that  taking  q  >  p  +  1  will  not  increase  the  robustness  rate. ) 

Proof  of  Theorem  4.1.  We  illustrate  the  proof  for  p  =  2,  for  the  two 
cases  q  *  2,  q  £  3.  By  Lemma  4.1,  we  see  that  for  both  cases,  for  k  £  1, 

(4.4)  Fq(N)  =  0|-«in(2,k-l)/2 

and  for  k£3  (i.e.  k  £  p  +  1),  we  have  the  best  rate  that  we  can  expect, 

i.e.  0(N  *),  We  now  calculate  the  robustness  rate  actually  achieved,  given 
by  max(FQ(N),  g(N))  where  g(N)  is  defined  by  (3.8). 

For  q  =  2,  we  see  that  using  (4.3),  we  have 

(4.5)  g(N)  =  sup  inf  Ju  -z|_  • 

per  2,1 

By  Lemma  4.2,  for  any  k  £  0,  we  see  that 

Cx  s  g(N)  s  C 

and  hence,  for  k  £  1,  the  robustness  rate  is  max(FQ(N),  g(N))  »  C.  Hence, 
this  method  is  not  robust.  By  (4.4),  for  k  £  3,  FQ(N)  *  CN  1,  so  that  by 
Theorem  3.3,  the  locking  is  0(N*). 

For  q  £  3,  we  get  (4.5)  again,  except  z  is  now  in  By  Lemma 

4.2,  for  k  £  1, 
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Cjh  S  g(N)  =s  C2h. 

i.e.  g(N)  «  CN-1/2.  Using  (4.4),  we  see  that  for  k  a  2,  FQ(N)  s  CH~1/2, 

-1/2 

so  that  the  robustness  rate  is  maxiF^CN),  g(N))  *  CN  for  k  a  2.  Also, 

for  k  a  3,  Fq(N)  *  CN  *  and  by  Theorem  3.3,  for  k  a  3,  the  locking  is 
0(N1/2). 


Remark.  For  pas,  the  above  results  hold  even  for  the  case  that  the 
quasiuniform  version  of  is  used  (see  [5]). 

Let  us  now  consider  the  uniform  rectangular  mesh  M^.  Suppose  first  that 

the  extension  y  is  based  on  ?>2(S)  elements,  with  *  V2’*1  = 

2  P  P<<1 

p*2’!?  x  ?>2’„.  Then  it  may  be  observed  that 
l  P.OJ  <1.0 


(4.6) 


W*  = 

q.i 


q  s  p 


®1»1*  y>2,h 

P+1.1  P,1 


q  a  p  +  1. 


The  following  is  an  analog  of  Lemma  4.2  for  this  case.  The  non-periodic 
version  of  this  result  is  established  in  [8],  [5]. 


Lemma  4.3.  Let  be  the  uniform  rectangular  mesh  of  Figure  3.1, 

for  u  €  Hlc+2(Q) ,  k  a  0,  q  a  1 , 
per 


Then 


(4.7) 


inf  |»  -z|2  * 

zeP2’? 

q.i 


Moreover,  there  exists  a  function  Q  e  C  (£3)  satisfying 

pe  r 


(4.8) 


zeP1’!?  ,u?>2’^ 
q+1.1  q, 1 


|Q  -zj.  a  Chq_1. 


Obviously,  the  two  bounds  (4.7)  and  (4.8)  hold  for  both  the  spaces 
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Now  using  Theorem  3.3,  Lemma  4.1,  (4.6)  and  Lemma 


t  u  TZ,h  and 
q+1,1  q, 1  q, 1 

4.3,  we  obtain  the  following  theorem,  whose  proof  is  similar  to  that  of 
Theorem  4.1. 


Theorem  4.2.  Let  the  extension  procedure  f  consist  of  the  h  version  on 
the  uniform  rectangular  mesh  M^,  with  spaces  {V^>  =  ^p*qj’  * 

|[*5]  x  ft*}.  Let  the  solution  sets  be  Then  for  the  p  and  q 

— r 

shown  below,  SF  is  robust  with  uniform  order  N  when  k  a  2r  +  1  and 
shows  locking  of  order  when  k  a  p+1 . 


order  *  0(N  r) 


Locking  order 
*  0(N*) 


We  now  consider  an  extension  procedure  on  M*1  based  on  :P3(S)  type 

u  3  h 

elements.  Accordingly,  suppose  v  *  V^*  .  Then  it  can  be  shown  that 


(4.7) 


q  s  p 


P+1.1 


p  *  2,  q  a  p  +  1,  or  (p,q)  *  (2,3). 


y>1,h  <j  j>2,h 

*3,1  U  r2, 1 


p  =  2,  q  a  4. 


In  this  case,  we  have  the  following  lemma. 


Lemma  4.4.  Let  be  the  uniform  rectangular  mesh  of  Figure  3.1  Then 

foi  «  €  H**^(Q),  k  a  0,  q  a  l 
per 


(4.8) 


Inf  |u 
*  b 


«  -  m, »in(k,m(q)V  . 

Az  *  Ch  Mk+2  . 


where  m(q) 
satisfying 


max(0,q-3).  Moreover,  there  exists  a  function  Q  €  C'  (Q) 

per 


(4.9) 


inf  |q 

h 


-  A,  *  Ch" 


Proof.  (4.8)  follows  from  Theorems  1,  2  of  [8],  as  shown  In  15).  Next, 

we  note  that  In  Lemma  5.2  of  [5J,  (4.9)  was  established  for  instead  of 

d.  1 

2  jj 

P * j,  for  the  non-perlodlc  case.  Essentially  the  same  proof  can  be  used  to 
prove  (4.9)  here  as  well.  □ 

Note  that  the  bounds  in  (4.8),  (4.9)  will  hold  for  IP  ,  and  T  ., 

d.l  d.l 

respectively,  as  well,  Using  (4.7)  and  Lemma  4.4,  we  may  once  again  establish 
the  following  theorem,  analogously  to  Theorem  4.1. 


Theorem  4.3.  Let  the  extension  procedure  f  consist  of  the  h  version  on 
the  uniform  rectangular  mesh  M^,  with  spaces  *  {Vpq}  * 

{[*5]  x  Let  the  solution  sets  be  d-  Then  for  the  p  and  q 

shown  below,  f  is  robust  with  uniform  order  N  r  when  k  a  2r  +  !  and 
shows  locking  of  order  N*  when  k  t  p  +  1. 
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Degree 

P 


Degree 

q 


Robustness 

Locking  order 

order  =  0{N  r) 

z 

o 

II 

r 

l 

0 

1/2 

0 

1 

1/2 

1/2 

(p-3)/2 

3/2 

(p-2)/2 

1 

Remark.  Theorems  4.2  and  4.3  show  that  locking  cannot  be  avoided  when 

rectangular  elements  are  used,  no  matter  what  choices  of  p  and  q  are  made. 

2  3  2 

Both  P  and  P  elements  are  robust  only  when  p  m  2  (for  ,  we  only 
P  P  Z 

3 

get  robustness  If,  moreover,  q  a  4).  For  p  i  3,  P  elements  show  twice  to 

2 

three  times  the  locking  rate  as  P  elements,  depending  on  the  choice  of  q. 


5.  THE  p  AMD  h-p  VERSIONS 


Let  us  now  consider  a  p  version  extension  procedure  P,  with  (v")  - 
{Vp'qJ,  where  h  is  fc®Pt  fixed  and  p.q  -*  ®.  Also,  let  us  consider  an  h-p 

version  over  a  quasiuniform  family  of  meshes  where  both  h  and  p,q 

are  changed  for  accuracy.  Then  we  have  the  following  estimate  for  FQ 
(Theorem  4.2  of  (51). 


Lsajaa  SLi.  Let 


<A  .  {/;*} 


be  a  sequence  of  p  version  spaces  on  a 


fixed  mesh  M*1,  with  p.q  -»  «.  Then  (3.3)  and  (3.4)  (for  any  fixed  d^  >  0) 


are  satisfied  with 


(5.1)  FQ(N)  *  CN  ‘,,c* 

with  C  independent  of  N  but  depending  on  h,  k.  Moreover,  as  p,q  -»  «#, 


-(k-l)/2 
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(5.1)  also  holds  If  the  h-p  version  over  a  quasiuniform  family  of  meshes  (M*) 
Is  used.  In  this  case,  we  have  the  following  more  refined  estimate: 


(5.2) 


F0(N) 


_  min(r, k-1 )  -(k-1 ) 
Ch  r 


where  r  *  mln(p,q). 

Note  that  In  the  above,  (M^>  does  not  have  to  be  a  family  of  uniform 
meshes,  but  can  be  a  quasiuniform  family  of  meshes. 

We  now  show  that  there  Is  no  locking  when  the  p  version  is  used  for  our 
model  problem. 


Theorem  5. 1 .  Let  the  extension  procedure  f  consist  of  the  p  version 
using  a  mesh  consisting  of  triangles  or  parallelograms  (which  can  be 
arbitrary).  Then  with  solution  sets  k  a  1,  y  Is  free  of  locking  and 

Is  robust  with  uniform  order  N  ^  as  p,q  -»  eo. 


Proof.  Using  the  results  of  [11],  [10],  we  can  show  that  for  the  p 

version  with  continuous  triangular  or  parallelogram 

per 

straight-sided  elements,  for  u  e  H^+2(Q),  k  a  1,  (h  fixed). 


per 


Inf  |u  -  z)  s  Cr’kMk,2 
r,  1 


2 

Now  for  any  1  ■  1,2,3,  with  V*1  =  V1,h  *  IV1’!?!  x  P1’*  we  have  T1, 

P.q  L  P*°J  r* 

1/2 

for  r  *  min(p.q).  51nce  r  *  0(N  ),  we  have,  with  g(N)  defined  by 

(3.8), 


w" 


g(N)  -  CN  , 

so  that  the  method  Is  robust  with  order  max(g(N) ,FQ(N) )  *  CN  by 

Lemma  5.1.  The  theorem  follows,  using  Theorem  3.3.  a 
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Let  us  remark  that  in  the  above  proof,  it  is  observed  that  g(N)  is  of  a 


smaller  order  than  FQ(N).  This  implies  that  for  the  limiting  case  (where  we 
have  the  biharmonic  problem),  the  p  version  actually  shows  am  increase  in  the 
rate  of  convergence  (by  one  order  of  p),  rather  than  a  decrease  due  to 
locking.  The  reason  is  that  the  solution  to  the  biharmonic  is  of  higher 
regularity  than  that  of  the  plate  problem,  amd  the  asymptotic  rate  of 
convergence  of  the  p  version  only  depends  upon  the  regularity  of  the  solution. 

For  the  h-p  version,  we  may  show  the  following  theorem  for  triangular 
meshes. 


Theorem  5.2.  Let  the  extension  procedure  f  consist  of  the  h-p  version, 
using  quasiuniform  meshes  consisting  of  triangles.  Let 

pi  5,  qsp  +  1.  Then  with  solution  sets  H,,  k  t  1,  y  is  free  of 


V1,h  with 

p.q 


locking  and  is  robust  with  uniform  order  N 
.k-1  -(k-1). 

h  p  ) . 


Yd’ 

-(k-1 )/2 


(or,  more  precisely. 


Proof.  The  essential  idea  is  to  use  the  following  estimate  for  the  h-p 
’  *  ^  continuous  triangular  elements  for  w  e  H^Cn),  k  a  1 


version  with  C 


per 


per 


(5.2) 


k  -k 


inf  |u>  -  z|_  *  ClT'r 


l.h 
r,  1 


M 


k+2 


provided  r  a  5,  r  as  k  ♦  1.  An  analog  of  (5.2)  has  been  established  in  (6) 
for  the  case  of  elements  —  a  similar  technique,  combining  (5.1)  with  a 

scaling  argument,  works  for  the  case  as  well.  For  V1*  =  with 

q  *  p  +  1,  p  fe  5,  we  have  c  "  .  Hence,  using  (5.2)  and  Theorem  3.3, 

we  have 

_  k  -k  _-k/2 
g(N)  =  Ch  p  =  CN 

in  (3.8).  The  theorem  follows  easily.  a 
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We  see  from  the  above  that  when  the  h-p  version  Is  used,  the  separate 


robustness  rates  of  the  h  and  p  versions  are  combined.  Theorem  5.2  therefore 
combines  the  results  of  Theorems  4.1  and  5.1,  and  says  that  If  the  p  version 
Is  combined  with  a  locking-free  h  version,  then  the  resulting  h-p  extension 
procedure  is  also  free  from  locking. 
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